Elliptic flow in nuclear collisions at the Large Hadron Collider 
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We use perfect-fluid hydro-dynamical model to predict the elliptic flow coefficients in Pb + Pb 
collisions at the Large Hadron Collider (LHC). The initial state for the hydrodynamical calculation 
for central A + A collisions is obtained from the perturbative QCD + saturation (EKRT) model. 
The centrality dependence of the initial state is modeled by the optical Glauber model. We show 
that the baseline results obtained from the framework are in good agreement with the data from the 
Relativistic Heavy Ion Collider (RHIC), and show predictions for the pr spectra and elliptic flow of 
pions in Pb + Pb collisions at the LHC. Also mass and multiplicity effects are discussed. 

PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Ld, 47.75.+f 



I. INTRODUCTION 

Azimuthal anisotropics of hadron spectra are a good 
measure of collective behavior of the dense particle sys- 
tem formed in the ultrarelativistic heavy ion collisions 
The origin of such anisotropies is thought to be 
the rescattering among the particles, which are ini- 
tially produced in isotropic partonic interactions. These 
anisotropies can be quantified by the second Fourier co- 
efficient V2, the so called elliptic flow coefficient, of the 
azimuthal hadron distribution. In non-central Au + Au 
collisions at the Relativistic Heavy Ion Collider (RHIC) 
at BNL the observed quite large v-i is regarded as one 
of the strongest signals pointing towards the formation 
of thermalized strongly interacting partonic matter, the 
Quark-Gluon Plasma (QGP). 

A common way to interpret the measured elliptic flow 
data is through full equilibrium, perfect-fluid hydrody- 
namical models, see e.g. reviews Jp, H, H|, with the 
Cooper-Frye freeze-out mechanism [5[. Consistency be- 
tween these type of models and the measured data at 
RHIC is regarded as evidence for low viscosity [a] and 
fast thermalization [?J of the QGP. However, a complete 
description of the low-px spectra and elliptic flow coeffi- 
cients of pions simultaneously with e.g. those of protons 
is problematic in this simple approach. 

To improve the modeling, several different approaches 
have been developed. If the chemical freeze-out is taken 
to happen before the complete kinetic freeze-out in the 
Cooper-Frye type of decoupling 8, 9] , a good simultane- 
ous description of e.g. the transverse momentum spec- 
tra of pions and protons can be obtained QJj]. How- 
ever, elliptic flow tends to be overestimated in these 
models. At the moment, perhaps the best agreement 
with the data is obtained from the hybrid models, which 
combine hydrodynamical treatment of the QGP and 
hadron cascade simulation of the hadronic interactions 
[Til H2I Il3l . flij . It is clear that a proper understand- 
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ing of hadronic interactions and the freeze-out mecha- 
nism are important in interpreting the observed data at 
RHIC. Also hydrodynamical codes with shear viscosity 
corrections have been recently developed for non-central 
collisions [la, lla, [12, UM, [l9( . These calculations indicate 
also low values of viscosity in the expanding matter. 

The key input to the hydrodynamical models are the 
initial state and the equation of state (EoS) of the QCD 
matter. When the initial state is obtained by fitting the 
hydrodynamical model to the RHIC data, the observed 
total hadron multiplicities fix the initial total entropy 
very well, but the transverse profiles of the initial densi- 
ties are not as well constrained, see e.g. Rcf. [17]. One 
possibility for controlling this uncertainty is to use theo- 
retically predicted initial conditions. In this approach it 
is possible to predict the initial state for different collision 
energies and nuclei as well. This has been our strategy 
e.g. in Refs. [U M, HI H3 . 

In [2(| we used the initial state from the EKRT fi- 
nal state saturation model [13] and showed that the per- 
turbative QCD + saturation + hydrodynamics approach 
gives a good description of the pt spectra of pions and 
kaons, and hadronic multiplicities in central Au + Au 
collisions at RHIC. We also presented the predictions for 
the hadron pt spectra in central Pb + Pb collisions at 
the Large Hadron Collider (LHC). We have also shown 
predictions for the elliptic flow at the LHC [Hill- The 
aim of the present work is to expand our previous studies 
of the hadron spectra and elliptic flow at the LHC and to 
probe the uncertainties in the predictions. A similar ap- 
proach with a calculated initial state is used in Ref. [26j , 
where the initial state is obtained from the Color-Glass 
Condensate model [27J. Other LHC predictions for the 
elliptic flow using the perfect-fluid hydrodynamics can 
be found in Refs. [25|, HH, [2^| and in the framework of 
viscous hydrodynamics in Ref. [3(3]. 

In this study we adopt the simple approach and use the 
perfect-fluid hydrodynamics for the space-time evolution 
of the matter, assuming full kinetic and chemical equi- 
librium throughout the evolution, see Sec. II. In spite of 
its restrictions discussed above, the full-equilibrium hy- 
ft.rodynamical approach, once tested against the RHIC 
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data, provides a good framework for baseline predictions 
of the low-px pion spectra and elliptic flow coefficients 
at the LHC. As discussed in Sec. Ill below, the EKRT 
model gives us the initial state for the hydrodynamical 
evolution in central A + A collisions for both RHIC and 
the LHC energies. Extension to non-central collisions is 
made by using the optical Glauber model at two differ- 
ent limits. The hydrodynamical evolution, elliptic flow, 
eccentricities and transverse flow at RHIC and LHC are 
discussed in Sec. IV. We will first, in Sec. V, show that 
the model results agree with the RHIC data, and then in 
Sec. VI present our predictions for the transverse momen- 
tum spectra of pions and for the elliptic flow coefficients 
of pions and protons at the LHC. Conclusions are given 
in Sec. VII. 



II. HYDRODYNAMICAL FRAMEWORK 

Once the initial energy density e and the net-baryon 
density ub are given and the Equation of State (EoS) 
P = P(e,n_e) is known, the evolution of matter can 
be described by relativistic hydrodynamics. By solving 
the perfect-fluid hydrodynamic equations, i.e. the local 
conservation of 4-momentum, d^T^ — 0, and of net- 
baryon number, ^(nsit* 4 ) = 0, one obtains the space- 
time evolution of all thermodynamic quantities and the 
collective flow velocity w' 1 = 7(1,vt,w z ) of the mat- 
ter. When considering particle production at midrapidi- 
ties, where the rapidity spectra are approximately flat, 
these equations can be simplified by assuming longitudi- 
nal boost invariance. In this case the longitudinal flow 
velocity is given by v z — z/t 31!], and all hydrodynamical 
quantities become independent of the space-time rapid- 
ity r\ — (1/2) In [(t + z) j (t — z)], i.e. they depend on 
the transverse coordinates x and y, and the longitudinal 
proper time r = \/t 2 — z 2 only. This reduces the (3+1)- 
dimensional problem to a (2+l)-dimensional one. We 
write then the conservation laws in these variables and, 
including the EoS, solve them numerically by applying 
the SHASTA algorithm [H, p| . 

In constructing the EoS [3J|, we describe the high- 
temperature phase as an ideal gas of massless quarks and 
gluons with number of flavors Nf = 3 and a bag constant 
B, while the low-temperature phase is taken as an ideal 
gas of all hadronic states with m < 2 GeV [35|. The QGP 
and hadron resonance gas (HRG) phases are connected 
via the Maxwell construction with mixed phase (MP) 
between the QGP and the HRG phase. The order of the 
phase transition has been shown to have only small effects 
on the hadron pr spectra and the elliptic flow of pions 
[36[ . As in our earlier studies [13, [Hj] , B is chosen such 
that the phase transition temperature is T c — 165 MeV. 
We assume full kinetic and chemical equilibrium for both 
phases throughout the temperature range considered in 
this work. 

Final hadron spectra are calculated through the 
Cooper-Frye decoupling procedure [5j as particle emis- 



sion from a constant-T surface obtained from the hy- 
drodynamic calculation. We determine the decoupling 
temperature Td ec ~ a parameter which controls the pt 
slopes of the hadron spectra in the single-Td ec hydrody- 
namic framework - from the RHIC data, as described in 
Refs. [13, [Hj]. For a recent discussion of the relation 
of Tdoc and decoupling dynamics, see [13, [H| . After the 
Cooper-Frye decoupling, all 2- and 3-body strong and 
electromagnetic decays of unstable hadronic states are 
accounted for. Thus, the feed-down from weak decays is 
not included in this work. 



III. INITIAL STATE AND CENTRALITY 
SELECTION 

As initial conditions, the boost-invariant hydrody- 
namic calculation, at an impact parameter &, requires 
the densities e(x, y, tq; b), ub{x, y, tq; b) at an initial time 
To. Our reference baseline is central collisions, for which 
we obtain the initial densities from the EKRT minijet 
(final state) saturation model [12]. In addition to the 
primary partonic transverse energy and the net-baryon 
number produced at midrapidity in central AA collisions, 
the EKRT model also gives the average formation time 
in terms of the saturation momentum, Tf = l/p sa t- 

Assuming immediate thermalization at production, 
To = Tf, using the binary collision (BC) profiles for 
the transverse-coordinate dependence of the densities, 
and setting Td cc = 150 MeV, we have previously shown 
that the minijet + saturation + hydrodynamics model 
is in good agreement with the RHIC data in most cen- 
tral Au+Au collisions [13, Hl|. We have also demon- 
strated that using the wounded nuclcon (WN) profiles 
and Tdec = 140 MeV leads to practically equally good 
results [37j . Also predictions for the central Pb+Pb col- 
lisions at the LHC have been presented [20I ]. 

The determination of the transverse profiles for e and 
riB is, however, problematic for the following reasons: 
First, the simplest version of the EKRT model for head- 
on collisions has only one saturation momentum scale, 
p sa t, which is perturbative, ^1 — 2 GeV for large nuclei 
at RHIC and LHC. In the localized version [§3] Ps&t de- 
pends on the transverse location but near the edges of the 
system, where the produced matter density becomes low 
enough, p sat becomes non-perturbative and the minijet 
calculation unreliable. Second, non-central collisions can 
be expected to constitute a multiscale problem as very 
different nuclear regions are colliding with each other, 
and we cannot expect the simple EKRT model to de- 
scribe the more peripheral collisions very well. Third, 
as discussed in [39|, [4(|, different regions in the trans- 
verse plane obviously form at different times; the center 
with a larger saturation scale forms earlier than the edges 
of the system. This phenomenon adds the complication 
of needing to provide the initial conditions on a surface 
t = tq(x, y) instead of just at a fixed initial time t = tq 
(see the discussion in |41|). 
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It is then obvious that some additional modeling is 
necessary in fixing the initial transverse profiles of e and 
riB for non-central collisions in particular. One possibil- 
ity is to assume a mixture of the BC and WN profiles, 
either for the energy density or the entropy density [4l| . 
with the mixture coefficient determined from RHIC data. 
Then, however, since the initial transverse profiles should 
be affected by production dynamics as discussed in [io| . 
an uncertainty in the extrapolation to the LHC energies 
would remain anyway since the mixture coefficients can 
change from RHIC to the LHC. 

Instead of introducing further model details and thus 
also further model uncertainties regarding the transverse 
profiles and their s/s dependence, we choose a simpler 
and more transparent approach: We fix the initial time 
To to the value computed from the EKRT model in cen- 
tral collisions, and we compute the spectra and elliptic 
flow coefficients with the BC and WN initializations sep- 
arately, considering these as two limiting cases of trans- 
verse profiles both at RHIC and at the LHC. 

Our BC initial conditions for A + A collisions are anal- 
ogous to the eBC model of Ref. [4l[ and are given by 



e B cO; b) 



AA 
BC 



(r;b) 



RC 



(r;b) 



1 [dE/dr,] b=Q dN } 
r N AA (0) dH 
1 [dN B /dr)] b=0 dNj A (b) 
r N4£(p) d 2 r 



(r;b), (1) 



where [dE/di]} b=Q and [dNB/dr]} b=0 and To are obtained 
from the EKRT model [2fj| for central collisions and 
where the BC profile is 



mil 

d 2 r 



r;b)=T A (r+~)T A (r-±)a% N , (2) 



and Ngg(0) 



a NN I ^ 2r pA(r)] 2 . In computing the 
standard nuclear thickness functions T A , we use the 
Woods-Saxon parametrization of the nuclear density. 
Similarly, our WN initial conditions are analogous to the 
eWN model of Ref. [4l[ and computed from 



e W N(r; b) 



(r;b) = d 



dN AA 
Q ttJV WN 

d 2 r 

dN AA 

flJV WN 



(Pi 



(r;b) 
(r;b) 



(3) 



where the normalization constants C e and Cb are fixed 
by requiring, for central collisions, the initial entropy 
dS/ drj and the initial net-baryon number in the eWN 
initial state to be the same as in the BC initial state. 
The WN profile is given by 



d 2 V 

= T A (r 



(r;b) 

+ 2 } 



+T A (r--) 



[1 
[1- 



A 



T A {v 



l-%T A (r + |) ]. (4) 



The initial state parameters for the different collision en- 
ergies are shown in Table I. The integral of dNy/N / d 2 r 
over the transverse plane is the number of participants, 
-^Vpart(b), for a collisions at impact parameter b. 



[GeV] 

to [fm] 

[dE/dr,] b=0 [GeV] 
[dN B /d V } b=0 
an K [mbl 



RHIC 

200 
0.17 
2460 
15.3 
42 



LHC 

5500 
0.097 
14800 
3.36 
60 



TABLE I: The initial state parameters for Au + Au collisions 
at RHIC and Pb + Pb at the LHC. 

In other words, we assume here that once the EKRT 
model fixes the normalizations of e and tib for central 
collisions, the optical Glauber model gives the impact 
parameter dependence of the initial conditions, assuming 
proportionality either to the BC or WN densities. We 
would like to emphasize that in this work we consider 
the two limits as they are. We do not try to fit the 
RHIC data by finding the best linear combination of the 
eBC and eWN profiles. Instead, we show that the RHIC 
data fall between the two limits and take the difference 
to represent the uncertainty in the extrapolation to the 
LHC energy. 

Qualitatively, there are two main differences between 
the two initializations. First, the number of BCs drops 
much faster with increasing impact parameter than the 
number of WNs. This leads to a faster dropping multi- 
plicity for the eBC than for the eWN initialization as a 
function of centrality. Second, at a given impact param- 
eter the BC density in the transverse plane falls faster as 
a function of the transverse distance from the center of 
the fireball than the WN density. At a given multiplic- 
ity, this leads to stronger transverse pressure gradients 
for the eBC than for the eWN initialization. Thus, if the 
same decoupling condition (same Tdcc) were used, the 
eBC initialization would lead to both stronger transverse 
flow and larger elliptic flow than the eWN initialization. 

To discuss the centrality classes, we apply the optical 
Glauber model, where the total cross section is given by 



A A 



= \ d 2 b 



dcr tot 
d 2 h 



= / d 2 b[l 



_ e T AA (h)a'£ N l 



(5) 



where T AA is the standard nuclear overlap function. For 
the centrality classes considered here, c\ = — 5%, 
c 2 = 5 — 10%, c 3 = 10 — 15%, we find the impact 
parameter ranges [0, 6i], [61,62], ... such that, e.g., for 
the 2nd centrality class c-i 



„AA /AA 
a c 2 /°tot 



0.05 



d 2 b[l 



T AA (b)<r- A 



(G) 



and similarly for the other classes c;. Using the dat t/d 2 b 
as the weight, we determine the average impact parame- 
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ter for each centrality class as follows 

(b)i = I ' d 2 b[l - e TAA ^ a ' l S»]b. (7) 

The average number of participants in a given centrality 
class is computed similarly. The initial densities for each 
centrality class are then computed using b = (6)j in Eqs. 
Q] and [31 The average impact parameters and average 
numbers of participants for selected centrality classes are 
shown in Table II. 



RHIC LHC 



centrality % 


b [fmj 




b [fmj 


A'part 


0-5 


2.24 


347 


2.31 


374 


5-10 


4.09 


289 


4.23 


315 


10-15 


5.30 


242 


5.47 


264 


15-20 


6.27 


202 


6.48 


221 


20-30 


7.49 


153 


7.74 


168 


30-40 


8.87 


102 


9.17 


112 


40-60 


10.6 


50.8 


10.9 


56.7 


60-70 


12.1 


19.6 


12.5 


21.2 


70-80 


13.0 


9.13 


13.4 


9.65 



TABLE II: The average impact parameters and the average 
numbers of participants in selected centrality classes from the 
optical Glauber model. 



pressure gradients stronger in the x direction. Stronger 
pressure gradients will generate stronger transverse flow 
during the evolution in the short direction. As a result 
the system grows faster in the x direction than in the 
y direction. The generated flow asymmetry manifests in 
the azimuthal dependence of hadron spectra and a non- 
zero V2 is observed. 

Figures 1 and 2 show the QGP-MP and MP-HRG 
phase boundaries and the decoupling boundaries in the 
x and y directions from our hydrodynamical simulations 
with b = 7.49 fm for RHIC and b = 7.74 fm for the LHC, 
both corresponding to the 20—30% centrality class. Both 
the RHIC and LHC results are shown with two different 
initializations, the eBC and eWN, discussed in the previ- 
ous section. Lifetimes of the different phases at different 
transverse locations can be read off from Figs. 1 and 2 in 
each case. With the same collision energy different ini- 
tial profiles correspond to slightly different multiplicity. 
Hence there is a difference e.g. in the initial transverse 
size of the eBC and eWN initializations in the figures. 
The first-order phase transition, present in the EoS we 
use, produces a shock wave at the boundary of HRG and 
MP, which in turn causes the peaks at the edge of the 
system seen especially in Fig. 2. However, contribution 
to the hadronic observables from the peaks is small. We 
have checked that removing the peaks causes less than 
5% changes in the hadron spectra and elliptic flow coef- 
ficients. 



IV. ELLIPTIC FLOW, ECCENTRICITIES AND 
TRANSVERSE FLOW 

The transverse momentum and rapidity dependent 
Fourier coefficients v n for each centrality class are de- 
fined as 

v n (y,p T ;b)=l—^-\ / ri0cos(n0) 
while the pr-integrated v n 's are given by 



(S) 



„„ (s;6)s ™- I r #c „ W) «, (9) 

V dV J J-7T dydcj) 

where b stands for the average impact parameter (b)i in 
the centrality class Cj. Due to the longitudinal boost 
symmetry assumed here, the v n coefficients, which we 
compute from the hydrodynamic spectra, do not depend 
on rapidity, thus y = is implicit. Our definition of the 
minimum bias elliptic flow coefficient, v 2 n ' hllis (y,PT), in 
turn is as follows, 



Jd%v 2 (y,p T ;b)j^ 



J clydpf^ 



(10) 



In the calculations the impact parameter is in the x di- 
rection. Therefore, the system is initially shorter and the 
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FIG. 1: (Color online) The phase boundaries and the 
constant-Td cc decoupling curves from our calculation for 
^/IWn = 200 GeV Au + Au collisions at RHIC in the 20- 
30 % centrality class. Both x and y directions are shown [l6| ]. 

A convenient way to illustrate the global features of 
the hydrodynamical space-time evolution of the matter 
is through the following three quantities: spatial eccen- 
tricity, momentum-space eccentricity and average trans- 
verse flow. These quantities are introduced in Ref. 0], 
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FIG. 2: (Color online) The Phase boundaries and the 
constant-Td e c MeV decoupling curves for -y/sjvjv = 5500 GeV 
Pb + Pb collisions at the LHC in the 20-30 % centrality class. 



where also their behavior for different systems has been 
discussed. The spatial eccentricity is defined as 



{y 2 - x 2 ) _ J dxdye(x,y,T)(y 2 - x 2 ) 
(y 2 + x 2 ) J dxdy e(x,y,T)(y 2 + x 2 ) 



(11) 



where the integral is over the transverse plane and the 
energy density e is used as the weighting factor. Similarly, 
the momentum space eccentricity is defined as 



/ dxdy (t xx - ryy) 

J dxdy (T xx + Tyy)' 



(12) 



where T M are the components of the energy-momentum 
tensor. The average transverse flow can be defined as 

_ (7U T ) _ / dxdye(x,y,T)jv T (x,y,t) 
[vt) = , ; = n—, — 7 s > 



(7) 



Jdxdye(x,y,T)j 



where vt — y v 2 + v 2 and the 7 factor has been used as 

an additional weight. 

The spatial eccentricity e x measures the asymmetry of 
the matter distribution in the x and y directions, thus it 
also measures the asymmetry of pressure gradients. In a 
noncentral heavy ion collision, with the impact parame- 
ter along the a;-axis, the initial distribution of the matter 
corresponds to a positive initial e x which acts as a driving 
force to an asymmetric transverse flow. The generated 
flow asymmetry is quantified by e p . The sign convention 
of e x and e p is such that an initially positive e x leads to 
a positive e p . A positive e p at the end of the evolution 
will convert to a positive elliptic flow coefficient V2 , which 



can be extracted from hadron spectra. The value of e p , 
when the system decouples, is approximately twice the 
PT-mtegrated V2 for pions [7(, see Figs. 3 and 5. During 
the evolution the established asymmetric flow field tends 
to drive e x towards zero, i.e. towards azimuthally sym- 
metric matter distribution. Therefore the driving force 
for the growth of elliptic flow is strongest during the early 
stages of the collision. At later times e x goes through zero 
and eventually to slightly negative values and e p will not 
increase anymore. 



Figure 3 shows the time evolution of e x , e p and 



from the same calculation as the curves in Figs. 1 and 2. 
Initially e x is very similar for RHIC and the LHC if the 
same initialization profiles (eWN or eBC) are used. At 
the beginning of the evolution, when most of the matter 
is in the QGP phase, the flow generated by the pressure 
gradients starts to decrease e x . Simultaneously e p is in- 
creasing. The initial rates of change for e x and e p are 
very similar at RHIC and the LHC. During the QGP 
phase there is no significant difference in the behavior of 
e x and e p between the two collision energies, but much 
bigger difference between the two chosen initial profiles. 
When the system enters the mixed phase all pressure 
gradients vanish ceasing the increase of e p and (vt), i.e. 
e p saturates even before e x goes to zero. During the 
mixed phase the matter is not accelerating, but the es- 
tablished flow field still expands the system transversally. 
This will cause s p actually to decrease during the mixed 
phase. The transition times between the different phases 
for different initializations can be read off from Figs. 1 
and 2. These times coincide with the structures seen in 
the behavior of e p , and (vt)- Transitions between the 
phases happen at different times at different transverse 
locations, therefore e.g. (vt) never completely saturates 
when most of the matter is in the mixed phase: part of 
the matter is always either in the HRG or in QGP phase, 
where the pressure gradients do not vanish. 

A significant difference between the LHC and RHIC is 
the lifetime of the QGP phase. At the LHC e p has more 
time to grow before the system enters the mixed phase, 
which leads to a larger e p at the end of the evolution. 
Therefore we also expect that the elliptic flow coefficient 
V2 will be larger at the LHC than at RHIC. Longer life- 
time of the QGP phase reflects also in the behavior of 
e x and (vt)- At RHIC e p saturates well before e x goes 
to zero, and e x is still positive when the system enters 
the HRG phase and e p is slightly increasing. Therefore 
there could still be some elliptic flow generated during 
the HRG phase at RHIC, i.e. t>2 would have some sensi- 
tivity on the decoupling condition. This can be verified 
by an explicit calculation and within our framework 



increases by ~ 20 — 30% at ATj 



part 



150 when the decou- 



pling temperature is changed from Td cc = 160 MeV to 
T dcc = 130 MeV. On the other hand at the LHC the sat- 
uration of e p and the sign change of e x happen more or 
less simultaneously and when the system enters the HRG 
phase e x is already negative and e p is slightly decreasing. 
At the LHC V2 is slightly decreasing when Td ec is de- 
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FIG. 3: (Color online) The spatial and momentum space 
eccentricity and the average transverse velocity as func- 
tion of the longitudinal proper time for ^/snn = 200 GeV 
Au+Au collisions at RHIC (thin solid and dashed curves) 
and ys^T = 5500 GeV Pb+Pb collisions at the LHC (thick 
solid and dashed curves) in the 20-30 % centrality class. 



creasing but the HRG effects are significantly smaller at 
the LHC than at RHIC. Thus the predictions for the pt~ 
integrated vi for the LHC are insensitive to the details of 
the HRG dynamics and therefore quite robust once the 
initial profile is fixed. 

In contrast to the behavior of e p , (vt) never completely 
saturates. Even at late stages of the collision (vt) can 



still increase significantly. Although the azimuthal asym- 
metry of the pressure gradients almost vanishes, the pres- 
sure gradient in the radial direction does not vanish. 
Even though the gradient is small in the HRG phase, 
additional transverse velocity is generated and can be 
easily seen in the slopes of hadron pr spectra. Thus the 
Pt spectra of hadrons remain sensitive to the decoupling 
temperature and HRG dynamics both at RHIC and the 
LHC. During the QGP phase the pressure gradients are 
strongest and the longer lifetime of the QGP phase at the 
LHC leads to a clearly larger transverse flow before the 
mixed phase than at RHIC. Regardless of the decoupling 
condition the pr spectra at the LHC are always flatter 
than at RHIC. 

The behavior of the pT-integrated i>2 for different ini- 
tializations can be quite easily determined from the be- 
havior of e p alone. The differential elliptic flow coefficient 
Vsipr) however depends not only on e p , but also on the 
slopes of the hadron px spectra, which in turn depend 
on (vt) and temperature. In general, increasing e p will 
increase t^fa r) a t fixed pr but increasing (vt) will de- 
crease it (42| ■ The net effect depends on the details of the 
HRG dynamics. Within our framework with full chemical 
and kinetic equilibrium in the HRG throughout the evo- 
lution, the net result is that i>2(pt) is quite insensitive to 
the decoupling temperature. At RHIC V2(pt) is slightly 
increasing and at the LHC slightly decreasing with de- 
creasing Tdcc- These changes are, however, small and 
practically within our framework, predictions for V2{pt) 
are independent of the decoupling temperature. Changes 
in the HRG dynamics would change this behavior, e.g . 
chemical freeze-out before kinetic freeze-out @, [fj [Io| 
would modify the dependence of spectral slopes on tem- 
perature and flow conditions and thus would affect the 
behavior of V2(pt) as a function of the decoupling con- 
dition. These effects are, however, not studied in this 
work. 



V. RESULTS FOR RHIC 

Figure 4 shows the calculated pion spectra for different 
centrality classes in ^/sjvjv = 200 GeV Au+Au collisions 
at RHIC, compared with the PHENIX data [13 . Shown 
in the figure are the results corresponding to the two dif- 
ferent initializations we consider here. The decoupling 
temperature is fixed to Tdcc = 150 MeV for the eBC ini- 
tialization, as in our previous works [20l . l2ll | . As seen in 
Fig. 3, the eWN initialization generates less transverse 
flow for a given decoupling condition, therefore - in or- 
der to reproduce the pion px spectra - the system must 
be allowed to decouple later than for the eBC initializa- 
tion. In the eWN case, T doc = 140 MeV describes the 
data well. With this difference in T dcc , both initializa- 
tions give an equally good agreement with the data for 
central and mid-central collisions. For more peripheral 
collisions the calculations start to separate and, as ex- 
pected, the eBC results fall below the eWN results but 
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FIG. 4: (Color online) The pr spectra of positive pions for 
-y/sjvjv = 200 GeV Au + Au collisions at RHIC compared with 
the PHENIX data The solid (dashed) lines are for the 

eBC (eWN) initialization. The centrality classes are indicated 
in the figure and are scaled by the increasing powers of 10~ . 
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FIG. 5: (Color online) The pT-integrated t>2 of charged 
hadrons for ^/sjvjv = 200 GeV Au + Au collisions at RHIC 
(thin lines) and ^/sjvjv = 5500 GeV Pb + Pb collisions at the 
LHC (thick lines) vs. the number of participants. The data 
is from the PHOBOS collaboration El. 



the data lies still well between the two limits considered. 
Fig. 4 also demonstrates how in most central collisions 
the hydrodynamically computed pion spectra reproduce 
the data well in the region px < 3 GeV, while the appli- 
cability region of hydrodynamic spectra in the 60-70 % 
centrality class is limited to pr < 1.5 GeV, only. 

Figure 5 shows our model calculations of the px- 
integrated v-i for charged hadrons as a function of the 
number of participants, compared with the PHOBOS 
data [44j |. For mid-peripheral collisions (iV part ~ 150 — 
200) the agreement with the data is quite good. How- 
ever, for the most central classes we underestimate the 
data - a typical feature also in several other hydrody- 
namical models. The underestimation of vi is usually 
associated with fluctuations in the initial geometry of 
the system (45j | , which are not accounted for in our treat- 
ment. For very peripheral collisions the model eventu- 
ally starts to overshoot the data but in these collisions, 
with the small system sizes and particle numbers, we 
expect, in any case, the validity of hydrodynamic mod- 
eling to deteriorate. From mid-peripheral to central col- 
lisions the eBC initialization generates more elliptic flow 
than the eWN initialization, as expected on the basis of 
Fig. 3. This results from the larger pressure gradient 
in the eBC initialization. However, for more peripheral 
collisions the behavior changes and the eBC results fall 
below the eWN ones. Since the total multiplicity for the 
eBC initial state drops faster with increasing impact pa- 
rameter, it also causes the lifetime of the QGP phase to 



drop faster. For very peripheral collisions the lifetime 
of the eBC system is considerably shorter than that of 
the eWN system. Thus the elliptic asymmetry will not 
have time to develop before freeze-out. Also the pres- 
sure gradients in the eBC case are larger than in the 
eWN case for central and mid-peripheral collisions, but 
become very similar for peripheral collisions. As a result, 
the elliptic flow predicted by the eBC initialization will 
eventually drop below the eWN results as we go towards 
more peripheral collisions. At RHIC the results coincide 
when 7V part ~ 70. In addition to the RHIC results, Fig. 5 
shows also the LHC prediction, which we shall comment 
on shortly. 

Figure 6 shows our results for the minimum bias vi(j>t) 
of pions in Au+Au collisions at RHIC (thin lines), com- 
pared with the RHIC data from STAR g| and PHENIX 
[4 71 ]. The agreement with the data is quite good up to 
Pt S 1-5 GeV for both initial profiles. As is the case 
with the pT-integrated «2, also V2{pr) is larger for the 
eBC than the eWN initialization from mid-peripheral to 
central collisions. However, for peripheral collisions also 
V2(pt) with the eBC initialization eventually falls below 
the eWN results for the reasons discussed above. This 
can be seen in Fig. 5 as the crossing of the eWN and 
eBC curves for RHIC. For the minimum bias V2(pt) all 
centrality classes between — 80% are included and the 
net result is that V2(pt) averages nearly to the same val- 
ues with both initializations at RHIC. 
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FIG. 6: (Color online) The minimum bias «2(pt) of pions for 
^/snn = 200 GeV Au + Au collisions at RHIC (thin lines) 
and ysjvlv = 5500 GeV Pb + Pb collisions at the LHC (thick 
lines) as a function of transverse momentum. The data is 
from the PHENIX and STAR [H collaborations. 
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FIG. 7: (Color online) Our predictions for the transverse mo- 
mentum spectra of positive pions for ^/sjvjv = 5500 GeV Pb 
+ Pb collisions at the LHC. The results with the eBC (eWN) 
initialization are shown by solid (dashed) lines. The centrality 
classes are the same as for RHIC in Fig. 4. 



VI. PREDICTIONS FOR THE LHC 

In Fig. 7 we show our prediction for the transverse 
momentum spectra of positive pions for different central- 
ity classes in ^snn — 5.5 TeV Pb+Pb collisions at the 
LHC. The charged hadron multiplicity in the most cen- 
tral collisions is ca. 2600 [2(j. Once the framework is 
tested in RHIC Au + Au collisions, the only freedom left 
is the choice of the initial profile and of the decoupling 
temperature Tdec- For the initial profile we consider the 
same limits as at RHIC, i.e. we do calculations with both 
the eWN and eBC profiles, with normalization fixed from 
the EKRT model for central collisions. As discussed in 
the context of dynamical decoupling [37j, we do not ex- 
pect Tdec to change significantly from RHIC to the LHC. 
Thus we show the eBC results with Td cc = 150 MeV and 
the eWN results with T doc = 140 MeV. Both initializa- 
tions give similar results for central and mid-peripheral 
collisions but, similarly to the RHIC case, the two calcu- 
lations start to separate for more peripheral collisions. 

In addition to the RHIC results, we show our predic- 
tions for the LHC in Fig. 5 for the pr-integrated v-i and 
in Fig. 6 for the minimum bias V2{pt), both for the eBC 
and eWN initializations. From Fig. 5 we see that, as ex- 
pected, the pT-integrated v-i is much larger at the LHC 
than at RHIC. Due to the longer lifetime of the QGP 
at the LHC the flow asymmetry has more time to de- 
velop before the system goes into the mixed phase where 
pressure gradients vanish. 



Although the pr-integrated t>2 is clearly larger at the 
LHC than at RHIC, the situation is not as clear for the 
minimum bias vi{j>t)- For the eBC initialization our 
LHC prediction is clearly above all the RHIC data points 
and, in the small-pr region, also above the RHIC predic- 
tions. On the other hand, in the low-p^ region the eWN 
results for the LHC are very close to the V2(pr) values 
at RHIC. This does not contradict the results for the pr~ 
integrated V2, which are obtained by weighting with the 
particle's px spectra. Since the spectra at the LHC are 
much flatter than at RHIC (,i.e. the average pr is larger 
at the LHC), i>2 gets more weight from the higher region 
at the LHC. Thus, even if the V2{pt) curves coincide, the 
PT-integrated V2 will be larger at the LHC. The difference 
in V2{pt) between the RHIC and the LHC results is that 
at the LHC the differential V2 is more sensitive to the 
initial profile than at RHIC, i.e. the eBC initialization 
leads to clearly larger values of V2(pt) than the eWN ini- 
tialization, while at RHIC both initializations give very 
similar results. The reason for this is the larger overall 
multiplicity and thus longer lifetime of the QGP phase at 
the LHC. Even in more peripheral collisions the lifetime 
of the QGP phase is still long enough to generate enough 
elliptic flow such that the ordering of t>2 (pt) between the 
eBC and eWN initialization remains down to N palt ~ 50. 
This can be seen in Fig. 5. 

Although our simple single- Td ec freeze-out model can- 
not reproduce in detail the measured proton pr-spectra 
and V2 at RHIC, it is still interesting to compare the gen- 
eral features in the calculated results for RHIC and the 
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FIG. 8: (Color online) The minimum bias V2(pr) of protons 
for ^sjvjv = 200 GeV Au + Au collisions at RHIC (thin lines) 
and v /s^zv = 5500 GeV Pb + Pb collisions at the LHC (thick 
lines) vs. p T . The RHIC data are from the STAR and 
PHENIX g7|, S| collaborations. 



LHC. Figure 8 shows the results from the same calcula- 
tion as for Fig. 6, but now for protons. Also the RHIC 
data from STAR g| and PHENIX HzL EH collabora- 
tions are shown. When the calculations are compared 
with each other, for both the eBC and eWN initial pro- 
files it is seen that the proton V2{pt) at fixed px at the 
LHC is always below the RHIC values. This can be un- 
derstood as a mass effect [Hj]. For the pions, the increase 
in average px compensates the increase in the flow asym- 
metry. However, for protons the increase of average px 
is larger due to the larger mass of the proton. Thus the 
increase of (pt) overcompensates the increase in the flow 
asymmetry. Although we cannot make a quantitative 
prediction of V2(pt) for the protons, we can infer from 
the results and arguments above that V2(pt) at the LHC 
is expected to be below the values measured at RHIC, 
while, as in the case of pions, the pT-integrated V2 is still 
expected to be larger at the LHC than at RHIC. 

Finally, we study the uncertainty from multiplicity in 
predicting the behavior of elliptic flow coefficients from 
RHIC to the LHC. Motivated by the multiplicity predic- 
tions in Refs. [IE 13, we have repeated the calculation 
for the LHC with half the multiplicity predicted by the 
EKRT model, by adjusting 



to 

[dE/d v ] b=0 



V2t 
1 

2V2 



[dE/d V ] b=0 , 



(14) 
(15) 



shows the integrated V2 for the two different multiplicities 
at the LHC together with results for RHIC and Fig. 10 
the minimum bias V2(pt)- With the lower multiplicity 
the lifetime of the QGP phase decreases and less elliptic 
flow is generated, which is clearly seen in Fig. 9, where 
the integrated V2 drops halfway between the RHIC re- 
sults and the original LHC prediction. The drop is sim- 
ilar between the eBC and eWN initializations. If a low 
multiplicity is observed at the LHC, the prediction with 
the eWN profile for V2 extends down to similar values 
as observed at RHIC. The integrated «2 is clearly sen- 
sitive to the observed multiplicity. The minimum bias 
V2(pt), shown in Fig. 10 is not as sensitive to the multi- 
plicity as the integrated elliptic flow, but there is still a 
visible difference as compared with the high-multiplicity 
prediction. Already for the high multiplicity calculation 
with the eWN profile, the LHC prediction is very close to 
the RHIC data and our calculations for RHIC. When the 
LHC multiplicity is lowered the eWN result drops at low 
Pt even slightly below the RHIC calculations and data. 
Also the eBC prediction drops at low pt closer to the 
RHIC results. 
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as suggested by the saturation conjecture [221 ]. Figure 9 
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FIG. 9: (Color online) The integrated V2 for the LHC low- 
multiplicity calculation 



VII. CONCLUSIONS 

We study elliptic flow in ultrarelativistic heavy ion col- 
lisions at RHIC and LHC energies using the framework 
of perfect-fluid hydrodynamics. The initial state for cen- 
tral heavy ion collisions is calculated from the pQCD 
+ saturation model [HI and perfect-fluid hydrodynam- 
ics is used to model the space-time evolution of the ini- 
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FIG. 10: (Color online) The minimum bias V2(pr) for the 
LHC low-multiplicity calculation 

tially formed matter. Centrality dependence is modeled 
by the optical Glauber model. We have shown that the 
results from this approach are in reasonably good agree- 
ment with the measured data at RHIC, and thus provide 
a good basis to predict both the pion spectra and elliptic 
flow coefficients in t/snn = 5.5 TeV Pb + Pb collisions 
at the LHC. The main uncertainty in our model is the 
initial transverse profile of the energy density. This un- 
certainty is addressed by considering two limits of the 
Glauber model, taking the energy density proportional 
either to the density of binary collisions or wounded nu- 
cleons. Three main observations are made for the LHC. 
First we note that the py-integrated vi is expected to be 
much larger at the LHC than at RHIC, since the lifetime 
of the QGP phase will be longer and thus the pressure 
gradients that generate elliptic flow will persist longer. 
The second observation is that «2(pt) for pions will be 
unchanged or may increase slightly from the values mea- 
sured at RHIC. This can be understood as an interplay 
between the flow asymmetry and flow magnitude [42j |. 
As the third point we observe that for heavier particles 
like protons, v-iipr) will be below the values measured 
at RHIC, even if the pT-integrated V2 is larger. This is 
due to the fact that the increased flow velocity is more 



important for protons, and thus it overcompensates the 
increased flow asymmetry. 

We further noted that the px integrated has still 
some sensitivity to the decoupling condition and there- 
fore to the HRG dynamics at RHIC, i.e. that spa- 
tial eccentricity is not entirely converted to the momen- 
tum space eccentricity before the matter enters the HRG 
phase. On the contrary, the prediction of v 2 at the LHC 
is more robust than at RHIC, once the initial profile 
and multiplicity are fixed: we do not expect the HRG 
dynamics to affect significantly the integrated v-i at the 
LHC. The differential elliptic flow coefficient v^pt) was 
in turn found to be insensitive to the decoupling con- 
dition at both collision energies. This insensitivity pre- 
sumably follows from our full equilibrium treatment of 
the HRG phase. Different HRG dynamics would change 
this behavior. However, when the decoupling condition 
is fixed at RHIC to reproduce the hadron spectra, the 
freeze-out temperature is not expected to change signif- 
icantly from RHIC to the LHC ,37?]. Once the EoS and 
the freeze-out parameters are fixed at RHIC, we expect 
that the same parameters apply also at the LHC. Once 
the multiplicities are fixed, the initial transverse profile 
remains as the largest uncertainty in the model. 

The effect of multiplicity on the elliptic flow coefficients 
at the LHC is studied by considering half the multiplic- 
ity that is predicted by the EKRT model. It is found 
that the integrated vi is quite sensitive to the multiplic- 
ity. The lower limit at the LHC, predicted by the eWN 
initialization, comes even slightly below the upper limit 
at RHIC given by the eBC initialization. The differential 
V2{pt) is found to be less sensitive, but the upper limit 
given by the eBC model, is moved closer to the RHIC 
results. 

We also note that we expect the applicability region 
of hydrodynamical modeling to extend higher in px at 
the LHC than at RHIC, because pions from jet fragmen- 
tation should start to dominate particle production over 
the hydrodynamical spectra at higher px at the LHC [2(| • 
Therefore we expect that the minimum bias i>2(pt) will 
reach higher values at the LHC even if the low-pr results 
are similar to those measured at RHIC. 
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